COUPLED PAINLEVE II SYSTEMS IN DIMENSION 4 AND THE 

SYSTEMS OF TYPE ^4^^) 
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Abstract. We find and study coupled Painleve II systems in dimension 4, which 
can be obtained by a degeneration from the systems of type ^4'-^''. We compare 
these systems with other types of coupled Painleve II systems from the viewpoint 
of the local index. We also give the phase spaces for these systems. 

0. Introduction 

In 2001, K. Takano described the confluence process in certain phase spaces of 
Painleve systems as deformations of manifolds [8]. In 1998, Noumi and Yamada [1] 

proposed the systems of type ^4^'' , which can be considered as a 4-parameter family 
of fourth-order coupled Painleve IV systems in dimension 4, whose Hamiltonians 
are given as follows: 

H = xy{2y - x - 2t) - 2l3iy - /^ax + zw{2w - z - 2t) - 2133W - {S^z + Ayzw 
= Hjv{x, y, t; P2) + Hiv{z, w, t; Ps, + 4:yzw. 

Here, x,y,z and w denote unknown complex variables, and ^1,^2,^3 and P4 are 
complex parameters. 

It is well-known that Pjv has a confluence to the second Painleve equation Pjj, 
where two accessible singularities come together into a single singularity. This sug- 
gests the possibility that there exists a procedure for searching for higher order 
versions of Painleve II, by using Takano's description of the confluence process from 
Pjv to Pji for the coordinate systems {x,y) and {z,w), respectively. In this vein, 
the goal of this work is to find a fourth-order version of the Painleve II equation. 
We present a 2-parameter family of fourth-order algebraic ordinary differential equa- 
tions which can be considered as coupled Painleve II systems in dimension 4, and 
which is given as follows: 



(1) 



dx 

Itt 
dy 

dt 

dz 

Itt 
dw 

~dt 



2 t 
-X + y + w - - 

2xy + ai 



-z +y + w-- 



2zw + a2. 
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Here x, y, z and w denote unknown complex variables, and a\ and a-i are complex 

parameters. 

Our differential system is equivalent to a Hamiltonian system. The Hamiltonian 
is given as follows: 



y2 t t 

H ^ — - {x"^ + -)y - aix + — - {z^ + -)w - a^z + yw 
= Eu{x, y, t\ ai) + Eu{z, w, t; + yw. 

From the viewpoint of symmetry, it is interesting to point out that our differential 
system is invariant under the following transformations: 
si : {x, y, z, w, t; ai, ^ (a; + ^, y, z, w, t; -ai, a2), 

52 : (x, y, z, w, t; ai, 0:2) {x, y,z + ^, w, t; ai, -a2), 

53 : {x,y,z,w,t;ai,a2) 

( rr -l+ai+a2+xy-yz l-ai-a2+xw-zw j.. -i ^,1 ^, '\ 
l^' x-z ^' ~z ' I, t - Q!2, i - ai), 

TT : {x, y, z, w, t; ai, 02) {z, w, x, y, t; 0:2, cti). 

We regard the system (1) as an algebraic vector field v on {x, y, z, w, t) e x B: 



d dx d dy d dz d dw d 
\ ^ ^ \ \ 

dt dt dx dt dy dt dz dt dw 



with B = C If we take a relative compactification P'^ x S of x B, the extended 
vector field v satisfies the following conditions: 

V e //°(p4,ep4(-logH)(H)). 

Here H is the boundary divisor in and <dpi(—logH){H) is the subsheaf of ©p4 
whose local section v satisfies f (/) e (/) for any local equation / of TC. Let us extend 
the regular vector field f on x 5 to a rational vector field v onF^x B. Then v has 
poles along the boundary divisor H. Moreover, v has accessible singularities along 
subvarieties in the boundary divisor 7i. (For the definition of accessible singularities, 
see Definition 1.1.) In order to explain our main results, we recall the definition of 
a symplectic transformation and its properties (cf. [4]). Let 

ip:x^xiX,Y,Z,W,t), y^yiX,Y,Z,W,t), z ^ z{X,Y, Z,W,t), 

w = w{X,Y,Z,W,t), t = t 

be a biholomorphic mapping from a domain D in 3 {X, Y, Z, W, t) into 3 
{x, y, z, w, t). We say that the mapping is symplectic if 

dx A dy + dz /\ dw = dX /\dY + dZ A dW, 

where t is considered as a constant or a parameter, namely, if, for t = to, (pt^ = <p\t=to 
is a symplectic mapping from the to-section Df^ of D to ip{Dt^). Suppose that the 
mapping is symplectic. Then any Hamiltonian system dx/dt — dH/dy,dy/dt — 

—dH/dx, dz/d,t = dH/dw, d,w/dt = —dH/dz is transformed to dX/dt — dK/dY, 
dY/dt = -dK/dX, dZ/dt = dK/dW, dW/di, = -dK/dZ, where 



dxAdy + dzAdw-dHAdt^dXAdY + dZAdW-dKA dt. 
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Here t is considered as a variable. By this equation, the function K is determined 
from H uniquely modulo functions of t. namely, modulo functions independent of 
X, Y, Z and W. For the vector field in (1), we can prove the following theorem. 

Theorem 0.1. The phase space X over B = C for the vector field v in (1) is 
obtained by gluing nine copies 0/ x C : 

C/o X C = X C 9 (x, I/, z, w, t), 
C/j X C = X C 9 {xj, yj, Zj, Wj, t) [j = 1, 2, 8), 

via the following symplectic transformations : 

1) xi = i, 1/1 = -x{xy + ai), zi = z, Wi = w, 

2) X2 = X, y2 = y, Z2 = W2 = -z{zw + a2), 

3) xs ^ x{l + {y + w - 2x^ - t)x - 2w{z - x) - ai - 0:2), Z/3 = ^, 

Z3 = x'^{z - ,x), = ^, 

4) xa = z^{x-z), y4 = ^, Z4 = z{l+{y+w-2z'^-t)z-2y{x-z)-ai-a2), W4 = i, 

5) X5 = 1/5 = -x{xy + ai), ^5 = i, W5 = + ^2), 

6) Xq = x{l + {y + w — 2x^ — t)x + 2{xw — zw — 0:2) — ai + 02), 1/6 = -, 

81 2)(a;'u;— zu;— 02) 



81 _ z"^ {z~x){yz-xy-ai) 



> 



2:7 = + (?/ + w - 2^;^ - + 2{yz - xy - ai) + ai - ^2), w-j = \, 
8) X8 = -((x - + ai + q;2 - l)y, ys = ^, = ^, ^8 = w + y. 

Because every coordinate transformation is symplectic, the Hamiltionian system 
in C/q X C is also written as a Hamiltonian system in each C/j x C (j = 1, 2, 7). 
By direct calculation, we can verify the following theorem. 

Theorem 0.2. On the affine open set {xi,yi, Zi,Wi,t) & Ui x B in Theorem 0.1, 
each Hamiltonian Hi on Ui x B is expressed as a polynomial in Xi,yi, Zi,Wi and t 
and satisfies the following conditions : 

dx A dy + dz A dw — dH Adt = dxi A dyi + dzi A dwi — dHi A dt. 

This paper is organized as follows. In Section 1, we consider the relation between 
the systems (1) and the systems of type ^4*^^^ In Section 2, we review the notion 
of accessible singularity and local index. In Section 3, we compare the systems (1) 
with other types of coupled Painleve II systems in dimension 4. In Section 4, we 
will prove Theorem 0.1 by giving an explicit birational transformation for each step. 
In the final section, we will prove Theorem 0.2. 



1. Relation between the systems (1) and the systems of type ^4^^^ 
As is well-known, the degeneration from P/y to Pn is given by 

1 ^ ^ -1 + e^T l + 2e^X eY 

Here the change of variables from (x, y) to (X, Y) is symplectic. As the fourth order 
analogue of the above confiuence process, we consider the following degeneration 
from the systems of type ^4^^^ to the systems (1). 
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The systems of type A^^^^ are given as follows: 



dx 

It 
dy 

dt 

dz 

It 
dw 

'dt 



-x^ + 4xy - 2tx - 2/3i + ^zw 
-2y^ + 2xy + 2ty + 
-z^ + ^zw - 2tz - 2/33 + %2; 
-2w^ + 2zw + 2tw + (34- Ayw. 



The systems of type ^4^''^^ are reduced to the systems (1) by putting 

X 1 1 I ^^^7"" 

A = TIT ' /52 = , /^s = TIT ' /54 = Q;2 , ^ = 



4e6 



V2e 



4e6 
^ V2 



1 + 2e2Z 



w 



V2e3 

71' 



and taking the limit e — > 0. 

2. Accessible singularity 

Let us review the notion of accessible singularity. Let S be a connected open 
domain in C and tt : W — > B a smooth proper holomorphic map. We assume that 
7i C W is a normal crossing divisor which is flat over B. Let us consider a rational 
vector field on W satisfying the condition 

V e //°(W,ew(-log7i)(7^)). 
Fixing to E B and P e Wt,,, we can take a local coordinate system (xi, ^2, x„) of 
yVto centered at P such that Ti smooth can be defined by the local equation xi = 0. 
Since v G if°(W, 6w(— log7i)(7i)), we can write down the vector field v near 
P= (0,0, ...,0,to) as follows: 



d d a2 d 



+ 



Xi dXr. 



This vector field defines the following system of differential equations 



' dxi 
"dt 
dx2 
"dt 



(2) 



ai(xi, X2, x„, t) 

a2{xi,X2, ....,Xn,t) 
Xi 



dXji (lfi{x\^ X2i — , Xfi^ t^ 

. dt Xi 

Here ai{xi, X2, t) {i — 1, 2, n) are holomorphic functions near P — (0, .., 0, to), 
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Definition 2.1. With the above notation, assume that the rational vector field v on 
W satisfies the condition 

V e 0w(-logH)(H)). 
We say that v has an accessible singularity at P — (0, 0, 0, to) if 

Xi = and 0,(0, 0, 0, to) = for every i, 2 < i < n. 

If P G Ttsmooth is not an accessible singularity, all solutions of the ordinary differ- 
ential equation passing through P are vertical solutions, that is, the solutions are 
contained in the fiber Wto over t = to. If P G 'Hsmootii is an accessible singularity, 
there may be a solution of (2) which passes through P and goes into the interior 

Let us recall the notion of local index. When we construct the phase spaces of 
the higher order Painleve equations, an object that we call the local index is the key 
for determining when we need to make a blowing-up of an accessible singularity or 
a blowing-down to a minimal phase space. In the case of equations of higher order 
with favorable properties, for example the systems of type A^-^^ [1], the local index 
at the accessible singular point corresponds to the set of orders that appears in the 
free parameters of formal solutions passing through that point [9] . 

Definition 2.2. Let v be an algebraic vector field which is given by (2) and {X, Y, Z, W) 
be a boundary coordinate system in a neighborhood of an accessible singularity P — 
(0, 0, 0, 0, t). Assume that the system is written as 



' dX 
"dt 


^a + fiiX,Y,Z,W,t) 


dY 


_ bY + f2{X,Y,Z,W,t) 


'dt 


X 


dZ 


cZ + fs(X,Y,Z,W,t) 


~dt 


X 


dW 


_ dW + f^{X,Y,Z,W,t) 


. dt 


X 



near the accessible singularity P, where a, b, c and d are nonzero constants. We 
say that the vector field v has the local index (a, b, c, d) at P if fi{X, Y, Z, W, t) is a 

polynomial which vanishes at P = (0, 0, 0, 0, t) and fi{X, Y, Z, W, t) {i = 2, 3, 4) are 
polynomials of order 2 m X, Y, Z, W. Here fi e C[X, Y, Z, W, t] [i = 1, 2, 3, 4). 

Remark 2.1. Wc arc interested in the case with local index (1, -,-,-) G Z"^. If each 
component of (1, ^, ^, ^) has the same sign, we may resolve the accessible singularity 
by blowing-up finitely many times. But when different signs appear, we may need 
to both blow up and blow down. 

3. Compare the systems (1) with other types of coupled Painleve II 

SYSTEMS in dimension 4 



It is known that certain reduction of the Drinfeld-Sokolov hierarchy of type ' 
reduces to a fourth order differential system with affine Weyl group symmetry 
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W{C2^^^) [2]. This system is an autonomous ordinary differential equations for 5 
unknown functions fo, fi, f2,Ui,U2 involving complex parameters aQ,ai,a2 satisfy- 
ing tto + 2ai + a2 = —4, which are given as follows: 



r 

dT 

dfi 
dT 
dh 
dT 

dui 

Ht 

du2 

It 



-2uifo - Oio 

{Ui - U2)h - Oil 
2^2/2 - Oi2 

{Ui + U2)Ui 



/1-/2 



— {Ui + U2)U2 — 



h 

fo — fl 
h 



where h' — 0, /o + 2/i + f2 ~ 4T -\-2huiU2. Setting the variables x, y, z, w, t and the 
nonzero parameter h as 

-ui fo U2 f^ + nrr h -i 

2 ' ^ 8 2 ' 8 

we can obtain the following coupled Painleve II systems in dimension 4: 



(3) 



— = 2xy + ai 



dx 

lit 
dy 

dt 

dz 

lit 
dw 

m 



—X + v + Sw — 7 
4 



—z + 3y + w - 



t 



—— — 2zw + 0:2. 



Here the Hamiltonian of this system is given as follows: 

H = — - {x"^ + -)y - aix + — - (z^ + -)w - a2Z + 3yw 

= Hii{x, y, t; ai) + Hii{z, w, t; 02) + ^yw. 

From the viewpoint of the local index, there are the following differences between 

these coupled Painleve II systems. 
Notation 

.{X,Y,Z,W) = {^ -2,^,1,^). 



Systems 


Accessible singularities 


Type of local index 


(1) 


(X,Y,Z,W) = (0,0,1,0) 


(-4,-1,-3,+!) 




(X,Y,Z,W) = (0,0.-2.0) 


(-4,-l,+3,-5) 


(3) 


(X,Y,Z,W) = (0,0,2,0) 


(-4,-l,-5,+3) 




(X,Y,Z,W) = (0,0,-3,0) 


(-4,-l,+5,-7) 


Pii X Pn 


(X,Y,Z,W) = (0.0.0.0) 


(-4,-1,-1,-1) 




(X,Y,Z,W)=(0,0,-1,0) 


(-4,-l,+l,-3) 
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Remark 3.1. / don't know whether the accessible singular points {X,Y, Z,W) = 
(0, 0, 2, 0) and {X, Y, Z, W) = (0, 0, —3, 0) of the systems (3) can be resolved or not. 

4. Proof of Theorem 0.1. 

In this section, we will give an explicit resolution of accessible singularities of the 
systems (1) and will construct a family of phase spaces for the systems. In the 
case of Painleve equations, we can obtain their phase spaces by adding subspaces 
of codimension 1 to the original space. But in the case of fourth order differential 
equations, we need to add codimension 2 subvarieties to the original space in addition 
to codimension 1 subvarieties. 

4.1. Accessible singularities of the systems (1) 

Let P be an accessible singular point in the boundary divisor 7i and {X, Y, Z, W) 
a coordinate system centered at P, where {X — 0} C H. Rewriting the systems 
in the local coordinate system (X, Y, Z, W), the right hand side of each differential 
equation has poles along Ti. If we resolve the accessible singularity P and the right 
hand side of each differential equation becomes holomorphic in the coordinate system 
(X', Y', Z', W) e [/ = C^, then we can use Cauchy's existence and uniqueness 
theorem of solutions. In order to consider a family of phase spaces for the system 
(1), let us take the compactification x i? of x B. Moreover, we denote the 
boundary divisor in P'^ by Ti. Fixing the parameter ctj, consider the product P'^ x 5 
and extend the regular vector field on C"^ x i? to a rational vector field v onF^ x B. 
The following Lemma 4.1 shows that this rational vector field v has six accessible 
singular points on the boundary divisor H x {t} C P^ x {t} for each t e B. 

Px P^ 



Figure 1 

Lemma 4.1. The rational vector field v has six accessible singular points Pi — 
{{X„Y„Z„W^) = (0,0,0,0)} {t = 1,2,3,4) and P^ = {{Xi,Yi, Zi,Wi)\Xi = Fi = 
Wi = 0, Zi = 1} and Pe = {(X3, F3, ^3, ^3)1X3 = F3 = Z3 = 0, W^3 = -1} ■ Here, 
(Xi,ri,Zi,l^i) = (i,f,f,f), {X,,Y2,Z,,W2) = (f,^if), {Xs,Y,,Zs,Ws) = 

(f'^f'f) (^4,l4,Z4,1^4) = (f,^^,^). 

This lemma can be proven by direct calculation. □ 

Remcirk 4.1. By the symmetry 

TT : {x, y, z, w; ai, a-i, aa) — > {z, w, x, y; ai, a^, a2), 




it is easy to see that 
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7r(Pi) = P2, = A- 



Now we are ready to prove Theorem 0.1. 

4.2. Resolution of the accessible singular point Pi 

In this subsection, we will give an explicit resolution process for the accessible 
singular point Pi by giving a convenient coordinate system at each step. 

By the following steps, we can resolve the accessible singular point Pi. 
Step 1: We blow up at the point Pi 

(1) V (1) (1) (1) ^1 

Ai Ai Ai 

Step 2: We blow up along the surface {{xi'^^\yi^^\ ^i*^^\ ■u;i*^^^)|xi'^^^ = yi'^^^ = 0} 



Xi 



0} 



Step 3: We blow up along the surface {(xi'^^), 1/1 2:1 '^^^ Wi*^^)) I l/i'^^^—ai = Xi^'^'^ — 



(2) 

xi(^)=xi(^), = .i^^^^.i^^ «;i(^)=^i(^). 

4.3. Resolution of the accessible singular point P5 

In this subsection, wc will give an explicit resolution process for the accessible 
singular point P5 by giving a convenient coordinate system at each step. 

By the following steps, we can resolve the accessible singular point P5. At first, 
we take the coordinate system {(X5, Y^,Z^,W^) = (Xi, Yi, Zi — 1, W^i)} centered at 

Step 1: We blow up along the curve {(X5, ^5,^5, ^5)1X5 = ^5 = W5 = 0} ^ 



^5 



(i)^W^ 

X5' 



Step 2: We blow down the 3-fold {ix5^^\y5^^\ Z5^^\w5^^^)\x5^^'^ = 0} ^ x 

.(1) 



Z5^^^ + 1 




Figure 2. W^-flop 
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The resolution process from Step 1 to Step 2 is well-known as P^-flop. In order to 
resolve the accessible singular point P5 and obtain a holomorphic coordinate system, 
we need to blow down the 3-fold Vi = x P^ along the P^-fiber. After we blow down 
the 3- fold Vi, we can resolve the accessible singular point P5 by only blowing-ups. 

Step 3: We blow up along the surfaces {{x5^'^\y5^'^\ Z5^'^\w5^'^^)\x5^'^^ — ys^^^ = 0} 
and {{x,^'\ys^'\z,^'\w,('^)\w,^'^ = z,('^ = 0} 

(2) (2) 
x,^'^^x,^^\ z,^'^-z,^^\ = 

Step 4: We blow up along the surfaces {{x5'^^\y5^^\ Z5'^^\w5^^'>)\x5^^^ = ys^^^ - 
a, = 0} and {{x,(^\ y^^'^ z,(^\ w,^'^)\w,(^^ - 0^2 = z,(^^ = 0} 

(4) (3) (4) 1/5^^^ - "1 (4) (3) (4) ""^S^^^ - «2 

4.4. Resolution of the accessible singular point Pq 

In this subsection, we will give an explicit resohition process for the accessible 
singular point Pq by giving a convenient coordinate system at each step. 

By the following steps, we can resolve the accessible singular point Pq. At first, 
we take the coordinate system {xs^'^\ys^^\ Zs~'^\w^^^^) centered at Pq. 
Step 1: We blow up at the point Pq 

(1) (1) (0) (1) ^8^°) (1) ^8^°) 

step 2: We blow up along the surface {{xs^^\ys''^\ z^^^\w^'^^^)\x^^^^ - z^^^^ = 
^8^^) = 0} 

^ (1) y (1) 

(2) •'^8 ~ ^8 (2) (1) (2) (1) (2) (1) 

Step 3: Weblowupalongthesurface{(x8^^\y8^^\-28^^\f^8^^^)|a^8^^^+Q;i+Q;2 — 1 = 

ys^''> = 0} 

XS^'^ ^ + X) " \ W^)=y8(^ .8(^)=.8(^ W^)=W^^ 

4.5. Resolution of the accessible singular point P3 

In this subsection, we will give an explicit resolution process for the accessible 
singular point P3 by giving a convenient coordinate system at each step. 

By the following steps, we can resolve the accessible singular point P3. 

Step 1: We blow up along the curve {(X3, F3, Z3, W^sjl^s = F3 = ^3 = 0} 

-^3 ^3 

Step 2: We blow up along the surface {{x3^^\y3^^\ Z3^^\w3^^^)\y3^^^ = xg^^^ = 0} 
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ys 



step 3: We blow up along the surface {{xa'^ ' ,y3^ > , > ,ws^ >)\x3'^ > = ws'^ > = 0} 



Step 4: We make a change of variables 



This change of variables is necessary for making the transition functions in the 
description of X symplectic [3]. It is easy to see that there are two accessible singular 
points Ps = {{xs'^'\y3^^\zs^^\ws^"^)\xs'^''^ +ws('^ = 2,1/3(4) = Ws'^'^ = 0,^3^') = 1} 
and Pe = {{xs^^\y3^'\ Zs('\ws^'^)\xs'^'^ + Ws^^^ = 2,ys^^^ = w^^^^ = 0,^3^"^ = -2} in 
the domain {ix3^^\y3^'^\ Z3^^\w3^'^^)\y3^^^ = 0} = C^. 

Step 5: We blow up along the surface {{x3^'^\ y3^^\ Z3^'^\w3^'^^)\y3^^^ — Z3^'^^ — 1 — 

0} 

X3^'^ = X3^'\ y3^'^ = y.^'\ z^^'^-"-^^, W^) = W^^ 

Step 6: We blow up along the surface {{x3^^\y3^^\ Z3^^\ W3^^'>)\x3^^'> -\- W3^^^ - 2 = 
^3^^) = 0} 

X3^'^ -^-^-^^^^^^^ , y3^'^-y3^'\ Z3^'^^Z3^'\ 

step 7: We blow up along the surface {{x3^^\y3^^\ Z3^^\w3^^^)\y3^^^ = -23 = 0} 

X3^'^ = X3^'\ y3^'^ = y.^'\ ^3^'^-^' = 
Step 8: We blow up along the surface {{x3^'^\y3^'^\ Z3^'^\w3^'^^)\x3^'^^ = ys^^^ = 0} 

X3(«) = ^, y^^'^-y^^'\ Z3^'^-Z3^'\ W«) = W^). 

Step 9: We blow up along the surface {(x3(^), 7/3^^), Z3^^\ W3''^^)\y3^^^ = Z3^^^ = 0} 

X3^'^-X3^'\ y3^'^-y3^'\ zs^'^-^, W^)=^3(«). 
Step 10: We blow up along the surface {{x3^^\y3^^\ Z3^^\w3^^^)\x3^^^ —t — y3^^^ — 

0} 



X3(^°) = ^^, y^^-ys^'K Z3^''^-Z3^'K = 

Step 11: We blow up along the surface {{x3^^^^\y3^^^^\ Z3^^"\w3^^^^)\y3^^^^ 
^3(10) _ 2^3(io)«;3(io) + 1 - «! - q;2 = 0} 
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4.6. Resolution of the accessible singular point Pq 

In this subsection, we will give an explicit resolution process for the accessible 
singular point Pg by giving a convenient coordinate system at each step. 

By the following steps, we can resolve the accessible singular point Pq. 

Step 1: We take the coordinate system centered at Pe 

Step 2: We blow down the 3-fold {{xG^^\y6^^\ Z6^^\ W6^^^)\y6^^^ = 0} C x x 

pi 

3 

Step 3: We make a change of variables 

Step 4: We blow up along the surface {{xQ^^\yG^^\ ZQ^^\we^^^)\xQ^^'> = ye^^^ = 0} 
Step 5: We blow down the 3-fold {{xQ^^\yQ^^\ ZQ^^\wQ^^^yQ^^^ - 0} C P^ x P^ x 

pi 

Step 6: We blow up along the surface {{xq^^\ yQ^^\ zq^^\ wq^^^)\x6^^^ = y^^^^ = 0} 

y6^^> 

Step 7: We blow down the 3-fold {ix6^^\y6^^\ Z6^^\w6^^'^)\ye^^^ = 0} C P^ x P^ x 

pi 

Xe^'^ = x,'-'K y,^'^=y,^'K z,('^ = y.'-^h,^^ , W^)=W'^. 

Step 8: We blow up along the surfaces {{xq^'^\ 2/6*•''^ -^6*•^^ W6^'^^)\x6^'^^ —t = y^^'^^ — 
0} and {{xe^'\ye^'\ze^'\we^'^)\ze^'^ = w^^'^ = 0} 

y6^^> ze(^> 

Step 9: We blow up along the surfaces {{xe'-^Kye'-^^ zq'-^\wq''^^)\xq'-^^ + 18?l'6^^^ + 
l-a^-a2 = i/e^^^ = 0} and {{xe^'\ ye^^\ ze^'\ we^^y)\we^^'> - f = ze^'^ = 0} 
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5. Hamiltonian system 
The system (1) is written as a Hamiltonian system: 

dx_dH^ dy__dH^ dz _ dH^ dw _ dH 
dt dy ^ dt dx ^ dt dw ^ dt dz ' 
Here the Hamiltonian is given as follows: 



H ^ — -{x + -)y - aix + — - [z + -)w - a2Z + yw. 



We list below the Hamiltonian for each coordinate system {xi,yi, Zi,Wi) {i — 
1,3,5,6). 

The Hamiltonian Hi in Ui We obtain 
1 



Hi — -(—twi + Wi + 2yi + txiyi — 2xiyiWi + Xiyi — 2ziWi + taiXi 



.2 , o.. , ^^-^^.^-^xlyiw ' -^-2 

—2aiXiWi + 2aixlyi + afxl — 2q;2-2i), 

where 



xi = - , yi = —xixy + cci) , zi^ z , wi^w. 

X 

The Hamiltonian Hi and coordinate system {xi,yi,zi,wi) above satisfy the fol- 
lowing condition: 

dxi A dyi + dzi A dwi — dHi A dt — dx A dy + dz A dw — dH A dt. 



The Hamiltonian H3 in U3 We obtain 

H3 = ^(2x3 - tys + yl + txsyl - 2x3yl + ^Ivt + '^ty^z^w^ - ^ylz^wj, + Axzylz2.W2. 

-2ylzlw3 + Aylzlwl + taiys - 2aiyl + '^(^i^^.vl + '^awlzi.wi, + alyl 
+to22/3 - 2a22/| + 202X31/1 - 2a22/|^3 + 4a2i/|^3tf^3 + 2aia2yl + alyl), 

where 

X XL) 

X3 — x{l+{y-\-w—2x'^—t)x—2w{z—x)—ai—a'2) , ys = — , z^ = x'^{z—x) , i/;3 = — . 

X X'^ 

The Hamiltonian H^ and coordinate system (0:3, 1/3, 2:3, 'u;3) above satisfy the fol- 
lowing condition: 

dxz A dys + dzz A dw^ — dH^ A dt — dx A dy + dz A dw — d{H + x) A dt. 



The Hamiltonian H^ in U5 We obtain 
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H5 = -{2w5 + 2^5 + txly5 + xlvl + tzlw^ + 2xly^zlw^ + z^wl + taix^ + 2aixlyr, 

^-2aix^zlw^ + a\xl + ^2^5 + 2Q;2a;52/52;5 + 2a2zlw^ + 2Q;iQ;2a;5^5 + «2^5), 
where 

3^5 = -, Db = -x[^y + (^i) ^ ^5 = -, W5 = + q;2). 

a: z 

The Hamiltonian i^s and coordinate system (xs, ^5, 2^5, Ws) above satisfy the fol- 
lowing condition: 

dx^ A dy^ + ^2:5 A dw^ — dH^ A dt — dx A dy + dz A dw — dH A dt. 



The Hamiltonian Hq in Uq We obtain 

He = ^(2x6 - tye + yl- l^^y^we + tx^yl - 2xeyl + ^Ivt - ^tyeZeWe 

+Aylz(iW6 - Axeylzewe + ^ylz^wl + taiye - 2aiyl + 2aiXeyl - Aaiylz^we 
+alyl - ta2y& + 2a2y6 - 2a2Xeyl + 4:a2ylzQWQ - 2aia2yl + alyl), 

where 

1 

Xq = x{l + + w — 2x^ — t)x + 2{xw — zw — Q!2) — cti + (^2) , ye = ~ i 

81 x^{x — z){xw — zw — a2) 

a;^(a; — z) 81 

The Hamiltonian Hq and coordinate system {xe,y6, ze,WQ) above satisfy the fol- 
lowing condition: 

dxQ A dye + dzQ A dwQ — dH^ A dt = dx A dy + dz A dw — d{H + x) A dt. 



The Hamiltonian ifg ^™ ^8 We obtain 

= ^{-tws + Wg - 2y8 + ixsyl - '^^IVs - + ^a^sZ/s^s - 2zlw8 + 2aiys 
-2aiXsyl + 2aiZ8 + ^a2y8 - ^0i2X8yl + 2a2Z8 - 2aia2y8 - 2aly8), 

where 

X8^ -{{x - z)y -l + ai + a2)y , = ~ ^ ^8 = ^, W8^w + y. 

The Hamiltonian iJg and coordinate system (xg, ys; -^S; ^s) above satisfy the fol- 
lowing condition: 
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dxs A dys + dzs A dws — dHs A dt — dx A dy + dz A dw — dH A dt. 
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